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INTRODUCTION 


A.l 

This addendum to the 50PACE 3-D (Boeing Plastic Analysis Canabllitv for 
3-Dimensional solids using isoparametric finite elements) orogram document 
(Ref.,, 1), describes modifications and additions incorporated into the orogram 
version dated September 1, 1975. These include the following major items. 

1. Restructuring of the program to allow solution of certain 3000-DCF 

v- 

problems within 64K words of core. (This also increases the 
allowable equation bandwidth/wavefront for larger core storage 
allocations). 

2. Generalization of the isoparametric element library, which now 
allows use of an element with any number of total nodes from 8 to 44, 
with any individual edge having from 2 to 5 nodes arbitrarily 
spaced along its length. 

i 

3. Generalization of the numerical integration scheme to allow an 
arbitrary number of Gauss -points (from 1 to 10) in each coordinate 
direction, with a default option which automatically selects the 
number of points in each direction such that the stiffness matrix 
will be exactly integrated. (This integration is generally exact 
only if the material properties are uniform within the element and 
all edges are straight lines). 

Reference 1. BOPACE 3-D, Final Report by The Boeinq Company to the NASA .v.,' 
Marshall Space Flight Center, Contract MASS-30615, April 15, 

1975. Document D180-13677-1 . 
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4. Provision of user-selected coordinate systens for describing the 
integration point locations and defining the integration point 
variables (stresses, strains, etc.). 


5. Elimination off the prohibition against disolacement constraints 
between nodes bn the same element. 

6. Addition of a capability for including elastic anisotropic 
materials, through user definition of the anisotropic thermal-strain 
behavior and stress-strain matrix, as functions of temperature. 


The following sections discuss the program modifications in more detail. Included 
are updates to the Drogran inout data formats, error messages, file usage, size 
limitations, and overlay schematic. /■ 
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A. 2 INPUT DATA 

> 

This section describes the items of ROPACE 3-D innut data which have been undated 
for the current version during incorporation of the new canabi lities. In oeneral , 
these updates are transparent to users of nrevious data decks, i.e. 1 , previous 
data decks can be run essentially without modification. 


C-3. The actual nroqram format for this item has Been corrected to aqree with 
the format given in the document. The correct format is (fiI5/3F10.Q) 
instead of (6I5,3F1O.0) . This now requires that inout item C-3 use 
two data cards instead of one (both cards may be blank if default values 
are used for the various constants in this item). 

1-1. Number of anisotropic materials (0-5), number of isotronic materials (0-5), 
fabrication temperature, and default values for four integration scheme 
codes. 

Format (6X, 212, VOX, F10.0, 30X, 14, 312) 

The number of anisotronic or isotropic materials may be zero, but the total 
number of materials must be from 1 to 10. 

The four integration scheme codes are defined as follows. 

1. The first code defines the tyoe of integration scheme, and in the 
current nrogram version may be given as either 0 (blank) or 1. 

Both of these current options result in product Gauss formula 
integration, with a variable number of points in each of the three 
parent £, n coordinate directions. The 0 specifies an auto- 

matic scheme, in which the program automatically selects the number 
of points in each direction such that the stiffness matrix will be 

original 

OF POOR QOM.KH a. 2-1 



exactly inteorated. (This intenration is qenerallv exact only if 
the material pronerties are uniform within the element and all 
edges are straioht lines)-, The 1 snecifies user-provided values 
for number of ooints in each direction, oiven as codes 2-4. 


2-4. 

pbiswab fam 
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These three codes are given if the first code has a value of 1. 
They are the number of noints (from 1 to 10) in each of the ^ _ 
5 , n»c coordinate directions, respectively. Total number 
of ooints may not exceed 600. It should be noted that the 
stiffness generation time for an element is highly deoendent 
upon the number of inteoration noints, as well as the number 
of element nodes. For examnle on the IBM 370, generation time 
for an 8-node element with 8 intenration ooints is less than 
1 second, while generation time for a 44-node element with 
125 integration points is almost 1 minute. 


1-2. This item now includes first the isotropic and then the anisotronic material 
data. (If the number of isotropic or anisotronic materials is 0, the 
corresnor.dinn data is omitted). The isotropic data is given as before. 

The anisotronic data is given •for each consecutive anisotronic material as 
described in the items below. 

a. Thermal strain curves. For each coordinate direction (in 

integration-point-variable coordinate system) give the following. 

Material number (inteoer from 6 to 10), coordinate component 
number (1, 2 or 3). 

Format (2110) 


Thermal strain curve for this coordinate direction. For each 
curve point give temperature and strain value, with points in orde^ 
of increasina temperature. User has option of from 1 to 4 noints 


ner card. 
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Format (8F10.0) 

Blank card after last point of each curve. 

b. Elasticity matrix definition. (Actual elasticity matrix relatinn 

total stress to total elastic strain, is computed at any temnerature 
as an Interpolated factor times the input matrix. The factor and 

input matrix are defined by the following). 

< / 

Material number (Integer from 6 to 10). 

Format (110) 

Curve of elastic’ matrix factors. For each curve ooint give 
temperature and factor, with Points in order of Increasing 
temperature. User has notion of from 1 to 4 points oep card. 

Format (8F10.0) 

Blank card after last point of curve. 

Elasticity matrix by rows, with six entries oer row. (Order is 
XX , YY, ZZ, XY, XZ, YZ). Matrix must be symmetric, and defines 
relationship between stresses and engineering strains. 

Format (6F10.0) 

For each element: element I.D. number, material number (1-5 for iso? 

tropic or 6-10 for anisotropic material), the eight corner nodes of the 
element, the maximum number of intermediate nodes along the edges, I.D. 
number of coordinate svstem to describe integration ooint locations 
( = 0, 1 or 2), I.D. number of coordinate system to define integration 
point variables ( = 0, 1, 2 or other), and the four integration scheme 
codes (if these codes are not given they are assigned the default values 
'from input item I -1 ) . 




Format (1015, 15, II, 14, 14, 312) 

If intermediate edge nodes are used, they are innut next for each edoe 
in the order shown in the figure (see Ref. 1). The same number of incuts 
per edge is expected (i.e., the maximum number of nodes ner edge as input 
on the previous card), with zeros used for edq&s which do not have the 
maximum number of intermediate nodes. Up to three intermediate nodes per 
edge may be used, and the total number of intermediate nodes may be un 
to 36. 

For each degree of freedom with a soecified displacement or constraint: 
give node I.D. number, component number (1, 2 or 3) and code. The code to 
be given is: 

1. For specified displacement, a blank (or may give the neaative of 

* 

the node I.D. number). 

2. For dependent constrained DOF, the node I.D. number of the 
independent DOF in constraint. 

User has option of specifying from one to four DOF on each card. 

Format (4(3I5,5X)) 

Blank card after last displacement-constraint DOF. 

Note: Specified forces may be denoted by giving the code eoual to the node 

I.D. number. However, this is not necessary because any DOF not soecified 
is automatically assumed to be a specified force. 



A. 3 Size Limitations 

The following variables are used to specify maximum size limitations in 
BOPACf 3-D. The values set for these variables are given in Table A. 3-1. The 
first set of values corresponds to the basic 64K proaram version, while the 
second set indicates the capabilities which would be obtainable throuah use 
of 128K core. 

NMAX1 - NMAX14 = values as defined in Ref. 1. 

MAXNPE = maximum number of nodes per element 

MAX I NT = maximum number of integration points per element 

B = maximum wavefront (in nodes) 
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Table A. 3-1: MAXIMUM PROGRAM SIZE LIMITATIONS 


VARIABLE 

64 K CORE 

128 K CORE 

Problem DOF 

3000 ® • 

3000 

NMAX1 

5 

5 

li ! ’ 

j 1000 

MMAX2 

1000 • 

NMAX3 

500 

500 

NMAX4 

5000 

5000 

NMAX5 

2000 

2000 

UMAX6 

20 

20 

i) 

NMAX7 

6 

6 

NMAX8A 

30 

30 

NMAX8B 

20 

20 

NMAX8C 

•c. 30 / 

30 

NMAX9 

10 

10 

NMAXTO 

6 

6 

NMAX1 1 

10 

10 

HMAXT2 

1000 

1000 

MM AX 13 

2 

2 

NMAX14 

60 

60 

MAXNPE 

44 

44 

MAX I NT 

600 

600 

B* 

61 

99 


* Based on MAIM program STQR array dimension of 38,000 for 64 K, and 102,000 
for 128 K. Larger STOP, array allows larger B. 



A. 4 ERRO& MESSAGES 

The following error message stoo code definitions have been added, deleted or 
modified. ' 

y>, •- 

35T- Wronq material number innut on an anistropic thermal strain card. 

352 Number of points input for an anisotropic thermal strain curve exceeds 
maximum. ’ 

355 No points input for an anisotropic thermal strain curve. 

// 

481 Wrong material number input on an anisotropic elastic matrix factor card. 

482 ^Number of points input for an anisotropic elastic matrix factor curve 

exceeds maximum. \ 

*>■ . . 

385 Mo points input for an anisotropic elastic matrix factor curve. 

718 I.D. of an integration point location coordinate system not equal to 0, 

1 or 2. 

719 I.D. of an intearation point displacement coordinate system is less th'an 
zero or greater than the maximum. 

721 The number of interior edae ; nodes on a side of the next element exceeds 3. 

720 The total number of integration points of an element is zero, negative, 
or exceeds the maximum. 

805 (Deleted) 
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FILE USAGE , 


A. 5 

The following files have been #i)ed to BOPACE. 


Unit Number 

Description 

Defined by 

INODAL (=21) 

Nodal data 

BLKDTA 

ISCP.l (=22) 

Contains temporary element data 
when BIGSC is called. Otherwise 

BLKPTA 


contains the load array PREF. 


ISCR2 (=23) 

Contains data that is to be added 
to the checkpoint tape 

\> 

BLKDTA 

i 



■A. 6 OVERLAY 

' 

The overlay of ROPACE was designed to minimize loading of seoments and to 
maximize the size of the MAIN orogram STOR array (common JLB) for a given 
core size. A schematic of the overlay is shown in Figure A. 6-1. 
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Figure A.6-1. Overlay Schematic 



